Graphene-based polaritonic crystal 
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It is shown that monolayer graphene deposited on a spatially-periodic gate behaves as a polaritonic crystal. 
Its band structure depending on the applied gate voltage is studied. The scattering of electromagnetic radiation 
from such a crystal is presented calculated and analyzed in terms of Fano-type resonances between the reflected 
continuum and plasmon-polariton modes forming narrow bands. 
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Coupling of light to the surface charges at a metal-dielectric 
interface gives rise to a special kind of evanescent electromag- 
netic (EM) waves called surface plasmon-polaritons (SPPs) 
|[T]. The specific properties of SPPs allow for their use in 
variety of practical applications. The sensitivity of SPPs to 
the properties of the dielectric, the metal and the interface is 
used in SPP-based sensors fT,^ and in high-resolution imag- 
ing mSl. Surface plasmons give rise to very large EM fields 
at the surface, which is important for surface-enhanced opti- 
cal spectroscopies |7|. Moreover, the SPP wavelength can be 
much smaller than the photon wavelength, opening the pos- 
sibility for further miniaturization of photonics components, 
a new field of research called nanoplasmonics f8|. Of par- 
ticular interest is the ability to tune SPP modes in plasmonic 
devices by external control: using an electric field in a liquid 
crystal |9|, a magnetic field in a magneto-optically active sub- 
strate 1 10|, thermal heating ifTTl . or a light beam focused on a 
non-linear coating lfT2l . 

The possibility of tuning the amount of free carriers in 
graphene using an external gate allows for an effective con- 
trol of the material's optical properties [[7311161 . Exploring 
graphene, a tunable two-dimensional (2D) metal, for plas- 
monics at the nano-scale reveals new physical effects and 
opens exciting possibiUties in this field lfT7l - [20l . When com- 
pared to their counterparts in conventional 2D electron sys- 
tems, SPPs in graphene exhibit some new and unusual prop- 
erties, such as the 1 /4-power density dependence of the SPP 
frequency fTT| and the existence of s-polaiized waves fTSl. 
Moreover, SPPs in graphene can potentially be used in a va- 
riety of practical applications. For example, using the ampli- 
fication of SPPs in graphene opens the possibility to create a 
terahertz radiation source ETl : employing the attenuated to- 
tal reflection (ATR) configuration with a gated graphene layer 
allows for a resonant switching of the reflection coefficient of 
an external EM wave from nearly unity to almost zero |22|. 

Recently, a new class of metamaterials based on SPPs in 
graphene, either deposited on a periodically corrugated sur- 
face 1 23 1 or composed of an array of micro-ribbons |24| was 
proposed, with a high potential interest for transformation op- 
tics. When SPPs propagate along a periodically modulated 
surface, the concept of a "surface polaritonic crystal" can be 
introduced [5^|, where the SPP dispersion shows a band-gap 
structure ll25l . in analogy with a photonic crystal. In this Let- 
ter we propose and theoretically analyze a new type of SPP 
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FIG. 1. (color on line) (a) Geometry of the system: "sandwich"-like 
structure containing a graphene layer on a substrate composed of a 
dielectric spacer (with dielectric constant ei) and a transparent con- 
ductive wafer with a periodic interface relief; (b) Real part of the SPP 
frequency versus wave vector in the first Brillouin zone, calculated 
for a single graphene layer at interface of two semi-infinite dielec- 
tric media; (c, d) Edges of the first two lowest gaps (at edge and at 
center of Brillouin zone) versus modulation depth a (panel c) or gate 
voltage V (panel d). The parameters are the following: £2 = 1, 
El = 3.9, /lo = 300 nm, D = 20.67 /^m, a = 0.6 (panels b and d). 
In panels (b) and (c) V — 50V, corresponding to an average chemical 
potential p. ~ 0.222 eV. In panel (b) the light lines for the cladding 
media, fcc/oj — Je\.2, are shown dashed (1) and dash-dotted (2). 



crystal, based on a graphene sheet deposited on top of a pe- 
riodically modulated gate electrode (wafer), as schematically 
represented in Fig. [TJa). 

In order to achieve the periodic modulation of graphene's 
conductivity, we consider a single graphene layer in the plane 
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z = deposited on a Si02 substrate with a dielectric con- 
stant £i [see Fig. [TJa)]- The opposite side of the substrate has 
a periodic relief with a spatial period D, such that h{x) — 
h{x + D). To be specific, we will consider the dielectric 
thickness modulation of the form h{x) = /io[l + acos{gx)], 
where a is the modulation depth and g ^ 2n/D. A con- 
ductive wafer is placed over this modulated surface, serv- 
ing as a gate contact. If a constant gate voltage V is ap- 
plied between the graphene layer and the wafer, a periodic 
modulation of the garphene conductivity can be achieved. If 
D 3> h{x), the carrier density in graphene can be expressed 
as n{x) = S^EiV / A'Keh{x), where e is the electron charge and 
^ is a coefficient between 1 and 2 depending on the charge dis- 
tribution on the corrugated gate electrode (further in the Letter 
we use ^ = 1). Owing to the periodicity of the surface relief, 
the carrier density in graphene is a periodic function with the 
same period D, thus resulting in the periodicity of its optical 
conductivity, a{x,Ljj) = a{x + D,uj). The latter can be re- 
lated to the local value of the chemical potential counted with 
respect to the Dirac point, fi{x) = hvp{nn{x)}^^^, where vp 
is the Fermi velocity. 

As it is shown below, the periodic modulation of the opti- 
cal conductivity leads to the possibility of direct coupling of 
a propagating EM wave to the surface plasmons (Fig. ^a), 
left, where £2 is the dielectric constant of the medium above 
graphene). However, only those SPP modes which lie within 
the light cone, u/k = c/y^, shown by dash-dotted lines in 
Fig. [TJb)' can be excited this way. In general, for the SPP 
excitation one has to consider an ATR structure like the one 
described in Ref. |22] which includes a prism with a dielectric 
constant £3 (Fig. [TJa), right). Usually there is a gap between 
the graphene sheet and the prism, which we shall model as 
a dielectric layer of a thickness d and a dielectric constant £2. 
We assume that the prism occupies the half-space z < ~d and 
a p-polarized EM wave impinges on the boundary z — —d, 
coming from z = —00 at an angle of incidence Q. 

Since the dielectric properties of the structure are periodic 
along x, the solution of Maxwell's equations, rotE*-'"-' = 
= — iKEmE*-™-' for the p-polarized wave 
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[with components E = {Ex,0, E^) and H = (0, Hy, 0)] can 
be written as Fourier- Floquet series: 

00 

H(^Hx,z)= J2 [Air^exp{nqi^'>z)+ (1) 

exp([«(fc + ng)x], 



n— — QQ 



si") exp(-Kq-'-z; 

00 (m) 

E(^^\x,z)^ J2 ^ exp(Kg(™)z)- (2) 

n— — oc 

B^") exp(-K(zi™)z)l exp([i(fc + ng)x]. 



In Eqs. 



and (2 1 



= \/[{k + ng)/^]^ - e,„ is the in- 
plane componenfof the photon wavevector in the medium m 
(to = 1, 2, 3; m = 3 applies only for the ATR configuration), 
K = Lo/c, and c is the speed of light in vacuum. Although 



the substrate is finite, for simplicity we shall consider the 
medium 1 as semi-infinite in the wave equations. If the gate 
electrode is transparent, this simplification is not crucial for 
the analysis of the optical properties of the modulated struc- 
ture. Boundary conditions at z — —d imply the continuity of 
the tangential components of the electric and magnetic fields, 
[Ei'\x,-d) - 
At z = 0, 
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the tangential component of the electric 
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field is continuous, Ei^\x,0) = Ex'^'{x,0), while the 
discontinuity of the tangential component of the mag- 
netic field, H^^\x,0) - H^^\x,0) = -(47r/c)> = 
-'{4:7r/c)a{x,uj)Ex{x,0), stems from the presence of 
surface currents (caused by the SPP electric field) in the 
graphene layer. Applying these boundary conditions, 
one can find the explicit form of the transfer matrices, 
j^pn^m+i^ which relate the coefficients in Eqs. ^ and 
§ for different to, (. . . , A^7\, 4+1, - 
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The matrices m"^™+i consist of 2x2 blocks M^;^™+\ 
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Snd is a Kronecker symbol, 

cr„ — Jf!^ a{x,uj) exp(—ingx)dx is the spatial 

Fourier transform of the graphene conductivity. 



The meaning of the coefficients An^^\ is different for 
the different media. The waves corresponding to the differ- 
ent terms in Eqs. ([T]l and (j2]i can be either propagating (with 
liii{q^^) < 0) or evanescent (with Re(g4'"^) > 0). Since the 
incident wave in the medium 3 possesses only the n = com- 
ponent (zero harmonic), Im(gQ ) < and the coefficients 
= Sn,oHiexp{—Kqj^^ d) are proportional to the mag- 
netic field amplitude (Hi) in the incident wave. In the medium 
1, the coefficients aI^^ = correspond to the absence of the 
corresponding harmonics coming from z = 00. Then, multi- 
plying the matrices Af ^^^M^'*"^, and taking into account the 
block-diagonal structure of the matrix AP^^, after some al- 
gebra we obtain the following equations for the amplitudes of 



the reflected haimonics, r„ = An^ exp{—Kqn'd): 
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FIG. 2. (color on line) (a) Zero harmonic reflection coefficient, 
\ro/Hi\^, vs frequency and incidence angle for SPPs in periodically 
modulated graphene with V = 50 V; (b-d) Zero harmonic reflection 
coefficient (red solid lines) and relative square amplitudes, \r„/Hi\'^, 
of n = 1 (blue dashed lines) and n = — 1 (black dashed lines) har- 
monics for the incidence angles S = 0° (b), O — 20° (c), and 
O — 80° (d) [the latter two correspond to the vertical dashed lines 
in panel (a)]; (e, f) Zero harmonic reflection coefficient vs frequency 
and gate voltage for 9 = 20° (e), and 6 = 80° (f). Other pa- 
rameters are the same as in Fig. [T| In panels (a, e, f) red (blue) 
color corresponds to high (low) values of the reflection coefficient. 
In panels (b, d) Fano-type fits are shown by dash-dotted (orange) 
fines, with Ro = 0.1508, SR = 0.035, q = 2, luq = lOmeV (b) 
and Ro = 0.116, SR = 0.08, q = 0.08, ujo = 10.68meV (d); 
7 = O.lmeV. 



where the elements of the matrix R and the vector HI, are 
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with Cn = cosh(Kg„ d) and Sn — sinh(K(7„ d). 

In order to obtain the general properties of SPPs in 
graphene with periodically modulated conductivity, we first 
consider the eigenvalue problem for the matrix R leading to 
the dispersion relation for SPPs in a flat 2D graphene layer 
placed between two lossless dielectric media (ei and £2)- We 
put £3 = £2 and solve the equation det(i?) = 0. It yields com- 
plex eigenvalues because the graphene conductivity a{x,Lu) 
has both real and imaginary parts, therefore the SPP eigen- 
modes are dissipative. The SPP dispersion curve for the real 
part of the frequency eigenvalue, uj', vs wavenumber for the 
first Brillouin zone, k £ [— g/2, g/2], is presented in Fig. 



[TJb). The imaginary part of the frequency (mode damping) 
is an order of magnitude smaller than a;'. As it can be seen 
from Fig. [TJb), the SPP dispersion curve is periodic in the 
fc-space, with the period g. There are bands of allowed SPP 
frequencies, separated by gaps opening at the edges and in the 
center of the Brillouin zone, where SPPs do not exist. As ex- 
pected, the widths of the gaps increase with the increase of 
the modulation depth a [Fig[TJc)]. A natural question arises: 
is it possible to control dynamically the gap widths through 
some nondestructive external knob? The positive answer to 
this question is evident from Fig. [TJd). Since the chemical po- 
tential of graphene can be tuned by the gate voltage, one can 
shift the spectral position and width of the gaps by changing 
V . Therefore, the SPP crystal band structure can be controlled 
dynamically. 

Another feature of the SPP spectrum in periodically mod- 
ulated graphene is that the "scan line", k — Ky^sinG , 
located within the dash-dotted lines in Fig. [ijb) crosses the 
SPP dispersion curves. This situation is completely different 
from the case of uniform graphene [22 1, where both phase and 
group velocities of SPPs are smaller than the velocities of light 
in the surrounding dielectrics. SPPs in periodically modulated 
graphene can be excited by an external propagating EM wave, 
without an ATR prism. This is illustrated by Fig. |2] where 
the amplitudes of the reflected field harmonics are presented; 
they have been calculated by solving Eq. (j3]l for d = and 
£3 = £2 (in this case, lo is real). 

At normal incidence (6 = 0), the zero harmonic reflection 
coefficient of the SPP crystal exhibits just one maximum at 
Lu « lOmeV [see Figs. |2|a) and (b)], which approximately 
corresponds to the upper edge of the second gap in Fig. [TJb). 
This is related to the parity of the SPP mode with respect to 
X ~ i), since it is excited by a plane wave and h{x) is an 
even function. Also from Fig. [2jb) it is clearly seen that the 
enhanced reflection of the zero harmonic corresponds to the 
excitation of the SPP harmonics with n = — 1 and n ~ 1. 
They correspond to the bottom of the third allowed SPP band 
and, for normal incidence, are mixed into the fc = band 
bottom mode. For oblique incidence [Figs. [2jc)-(f)], there 
are two resonances corresponding to the second and third SPP 
bands and producing reflected field harmonics with n = ±1. 
SPPs are effectively excited when the frequency and the in- 
plane component of the wave vector of the incident EM wave 
match those of SPP eigenmodes of the modulated graphene. 
The energy of the incident wave is transferred into the SPP 
harmonics with n = ±1, while the refleted wave in the far 
field contains only zero harmonic. Note that the position and 
the amplitude of the resonances can be controlled by changing 
the gate voltage [Figs. [2je) and[2jf)]. Quite interestingly, the 
direct excitation of n = ±2 and higher SPP bands can produce 
propagating EM waves with |fc ± < Ky^, scattered at 
angles 
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Let us focus on the characteristic spectral dispersion of the 
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FIG. 3. (color on line) Reflection coefficient of zero harmonic 
ro /-ffi P vs frequency oj and angle of incidence for ATR excitation 
of SPPs in graphene with parameters £3 = 16, £2 = 1, ei = 3.9, 
ho = 300 nm, D = 10.34 ^m, a = 0.6, d = 10 /im, V = 50V 
(panel a), or 90V (panel b). 



between the first and the second SPP bands. The mode an- 
ticrossing, corresponding to the edges of the gap, is clearly 
seen in Figs. |3ja) and[3];b) (near 8 « 50°, see insets). Com- 
parison of Figs. [3|a) and[3|b) shows, that increasing the gate 
voltage results in an increase of the gap width, as it could be 
anticipated from Fig. [TJf). 

To conclude, we have demonstrated that a single graphene 
layer deposited on a "sandwich" -like structure with a periodi- 
cally corrugated gate electrode has the properties of a polari- 
tonic crystal, namely, possesses a band structure with gaps 
that can be tuned by the gate voltage. This crystal exhibits 
Fano-type resonances in the reflectance of incident EM waves 
due to the excitation of plasmon-polariton Bragg modes. 



zero harmonic reflection coefficient around each of the SPP 
resonance frequencies [Figs. [2]^b) and (c)] where an intrin- 
sic mode of the SPP crystal is excited. This is an asymmet- 
ric Fano-type resonance with a maximum accompanied by a 
neighboring minimum, a phenomenon first discovered in the 
ionization spectra |26 | of He and now known in many differ- 
ent branches of physics including light scattering by photonic 
crystals |27 1. If a discrete mode couples to a continuum of ex- 
citations, the Fano resonance arises from the constructive (de- 
structive) interference of the localized and delocalized waves, 
taking place above (below) the discrete mode frequency uq. 
A simple analytic expression to describe such a lineshape was 
suggested by Fano ||26J : 



{qj + ui - uJo)'^ 



(6) 



where 7 is the mode damping and q is an assymetry param- 
eter related to the relative strengths of the transitions associ- 
ated with the discrete and continuum modes. Fig. ^h) shows 
a Fano-type fit, R{uj) = Rq + SR ■ F{uj), to the calculated 
reflectivity spectrum (i?o and SR are constants). The cou- 
pling between the continuum of propagating EM modes in the 
medium 2 and the SPP Bragg mode of the polaritonic crystal 
decreases with the increase of the incidence angle. When Q 
is close to the Brewster angle, 8^ = atan(-y/ei/e2) ~ 63°, 
the resonant excitation of n = ±1 modes results in two nar- 
row and almost symmetric minima in the spectrum [see Fig. 
[2|d) where the Fano fit corresponds to q — 0.08, compared to 
«7=.2inFig.|2];b)]. 

As mentioned above, the direct excitation of SPPs by a 
propagating wave allows for probing only the part of the po- 
laritonic crystal band structure comprised between the dash- 
dotted light lines in Fig. [ijb). Beyond this, one has to use 
the ATR scheme where SPPs are excited by an evanescent 
wave with a sufficiently large wavenumber, k = K^/e^ sin 8 . 
The results calculated for the ATR structure are presented in 
Fig. [3] where the SPP excitation conditions correspond to a 
minimum of the zero harmonic reflectance related to the gap 
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